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Abstract 

This paper explores integrable structures of a generalized melting 
crystal model that has two g-parameters qi,q2- This model, like the or- 
dinary one with a single g-parameter, is formulated as a model of ran- 
dom plane partitions (or, equivalently, random 3D Young diagrams). 
The Boltzmann weight contains an infinite number of external poten- 
tials that depend on the shape of the diagonal slice of plane partitions. 
The partition function is thereby a function of an infinite number of 
coupling constants ti,t2, ■ ■ ■ and an extra one Q. There is a compact 
expression of this partition function in the language of a 2D complex 
free fermion system, from which one can see the presence of a quan- 
tum torus algebra behind this model. The partition function turns out 
to be a tau function (times a simple factor) of two integrable struc- 
tures simultaneously. The first integrable structure is the bigraded 
Toda hierarchy, which determine the dependence on ii,t2,.... This 
integrable structure emerges when the g-parameters qi, q2 take special 
values. The second integrable structure is a g-difference analogue of 
the ID Toda equation. The partition function satisfies this g-difference 
equation with respect to Q. Unlike the bigraded Toda hierarchy, this 
integrable structure exists for any values of (?i , (?2 • 
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1 Introduction 



The melting crystal model is a model of statistical physics and describes a 
melting corner of a crystal that fills the first quadrant of the 3D Euclidean 
space. The complement of the crystal in the first quadrant may be thought 
of as a 3D analogue of Young diagrams. These 3D Young diagram can 
be represented by plane partitions. Thus the melting crystal model can be 
formulated as a model of random plane partitions. 

This model has been apphed to string theory P and gauge theory [21 
[3]. From the point of view of gauge theory, the partition function of the 
melting crystal model is a 5D analogue of the instanton sum of 4D Af = 2 
supersymmetric Yang-Mills theory [U El I^- (Curiously, the 4D instanton 
sum also resembles a generating function the Gromov-Witten invariants of 
the Riemann sphere [71 [8].) This analogy will need further explanation, 
because the 4D instanton sum is a sum over ordinary partitions rather than 
plane partitions. The fact is that one can use the idea of diagonal slicing 
to rewrite the partition function of the melting crystal model to a sum over 
ordinary partitions [2]. Comparing these two models of random partitions, 
one can consider the melting crystal model as a kind of q-deformation of the 
4D instanton sum. Here g is a parameter of the melting crystal model related 
to temperature. 

In our previous work [10] (see also the review [H]), we introduced a set 
of external potentials into this model, and identified an integrable structure 
that lies behind this partition function. Namely, the partition function, as 
a function of the coupling constants ti,t2, . . . of potentials, turns out to be 
equal to a tau function (times a simple factor) of the Toda hierarchy [121 IS] • 
Moreover, the tau function satisfy a set of constraints that reduces the full 
Toda hierarchy to the so called ID Toda hierarchy. Though a similar fact 
was known for the 4D instanton sum [131 [T5l [T6] , we found that the partition 
function of the melting crystal model can be treated in a more direct manner. 
We derived these results on the basis of a fermionic formula of the partition 
function [T3j. A technical clue is a set of algebraic relations among the basis 
of a quantum torus (or cylinder) algebra realized by fermions. These relations 
enabled us to rewrite the partition function to a tau function of the Toda 
hierarchy. 

In the present paper, we generalize these results to a melting crystal 
model with two g-parameters qi,q2 [17]. Actually, since the potentials have 
another g-parameter q, this model has altogether three g-parameters gi,g2 
and g; letting gi = g2 = g, we can recover the previous model. 

Our goal is two-fold. Firstly, we elucidate an integrable structure that 
emerges when gi and g2 satisfy the relations gi = g^/^i and g2 = g^/^^ fg]; 
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a pair of positive integers A^^i and The partition function in this case 
turns out to be, up to a simple factor, a tau function of (a variant of) the 
bigraded Toda hierarchy of type (A^i, N2) [18], which is also a reduction of the 
Toda hierarchy. Secondly, without such condition on the parameters qi,qi 
and q, we show that the partition function satisfies a g-difference analogue 
[ini [201 [211 [22] of the Toda equation with respect to yet another couphng 
constant Q. In the gauge theoretical interpretation, Q is related to the energy 
scale A of supersymmetric Yang-Mills theory. 

This paper is organized as follows. Section 2 is a review of combinatorial 
aspects of the usual melting crystal model. The model with two g-parameters 
is introduced in the end of this section. Section 3 is an overview of the 
fermionic formula of the partition function. After reviewing these basic facts, 
we present our results on integrable structures in Sections 4 and 5. Section 4 
deals with the bigraded Toda hierarchy, and Section 5 the g-difference Toda 
equation. We conclude this paper with Section 6. 

2 Melting crystal model 

In the following, we shall use a number of notions and results on partitions. 
Young diagrams and Schur functions. For details of those combinatorial tools, 
we refer the reader to Macdonald's book [23|- See also Bressoud's book [21] 
for related issues and historical backgrounds. 

2.1 Simplest model 

Let us start with a review of the ordinary melting crystal model with a single 
parameter g (0 < g < 1). As a model of statistical physics, this system can 
take various states with some probabilities, and these states are represented 
by plane partitions. 

Plane partitions are 2D analogues of ordinary (one-dimensional) parti- 
tions A = (Ai, A2, . . .), and denoted by 2D arrays 

... \ 

of nonnegative integers vTjj (called parts) such that only a finite number of 
parts are non-zero and the inequalities 



TT 



00 
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are satisfied. Let Ivrl denote the sum 



oo 



TT = 



i,j=0 



of all parts 7r.y. 

Such a plane partition vr represents a 3D Young diagram in the first 
quadrant x,y, z > of the (x, y, z) space. In this geometric interpretation, 
TTij is equal to the height of the stack of cubes over the (i,j)-th position of 
the base (x, ?/) plane. Therefore |7r| is equal to the volume of the 3D Young 
diagram. 

In the formulation of the melting crystal model, the complement of the 
3D Young diagram in the first quadrant embodies the shape of a partially 
melted crystal. We assume that such a crystal has energy proportional to 
|7r|. Consequently, the partition function of this system is given by the sum 



of the Boltzmann weight g''^' over all plane partitions vr. 
2.2 Diagonal slicing 

We can convert this model of random plane partitions to a model of random 
partitions by diagonal slicing. This idea originates in the work of Okounkov 
and Reshetikhin [9] on a model of stochastic process (Schur process). 

Let 7r(m) (m G Z) denotes the partition that represents the Young di- 
agram obtained by slicing the 3D Young diagram along the diagonal plane 
X — y = m 'm. the (x, y, z) space. In terms of the parts vTjj of the plane 
partition, these diagonal slices can be defined as 



The plane partition can be recovered from this sequence {7r(m)}^^_g^ of par- 
titions. To be diagonal slices of a plane partition, however, these partitions 
cannot be arbitrary, and have to satisfy the condition 




TT 




if m < 



m > 



< 7r(-2) -< 7r(-l) -< 7r(0) ^ 7r(l) ^ 7r(2) > . 



(1) 



Here 'V denotes the interlacing relation 



A = (Ai, A2, . . .) ^ = (yUl, /i2, • • •) 



Al > /il > A2 > /i2 > ■ ■ ■ , 
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namely, X >~ fi means that the skew diagram A//i is a horizontal strip. 
These diagonal slices, in turn, determine two Young tableaux 

T = {Tij)(ij)^x, T' = (Z;^)(ij)eA 
on the main diagonal slice 

A = 7r(0) 

as 

Tij = m if (i, j) G 7r(— m)/7r(— m — 1) 
T-j = m if {i,j) e 7r(m)/7r(m + 1). 

Since the skew diagrams '7r(±m)/7r(±(m + 1)) are horizontal strips, these 
Young tableaux turn out to be semi- standard tableaux, namely, satisfy the 
inequalities 

Tij > Tj+ij-, Tij > Tj (2) 

These semi-standard tableaux T, T' encode the left and right halves of the 
plane partition. The Boltzmann weight g''^' thereby factorizes as 

gi-i = gV, (3) 

where 



^ JJ ^(m+l/2)|7r(-m)/7r(-m-l)|^ ^T' = JJ g( 



_,(m+l/2)|7r(m)/7r(m+l) 

m=0 m=0 



Since the triple (A,T, T') is in one-to-one correspondence with the plane 
partition vr, the partition function Z can be reorganized into the sum over 
the partition A = 7r(0) and the sum over the pair (T, T') of semi-standard 
tableaux of shape A: 

z=i: E ^v. (4) 

A r,T':shapeA 

By the well known combinatorial definition of the Schur functions, the par- 
tial sum over semi-standard tableaux becomes a special value of the Schur 
function Sx{x) of infinite variables x = {xi, X2, ■ ■ ■) as 

E E ^'^' = ^a(/), (5) 

T:shapeA T'rshapeA 



^Note that the entries of the tableaux are arrayed in decreasing order rather than usual 
increasing order. This difference is immaterial as far as we consider Young tableaux on a 
fixed Young diagram A to describe the associated Schur function s\{x). 
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where 



Note that the g-weights q and q are identified with the monomials 



n 



(i,i)GA 



in the combinatorial definition 



Tishape A 

of Schur functions. We thus obtain the Schur function expansion 

Z = J2s,iqn' 
X 

of the partition function. 

2.3 Calculation of partition function 

We can calculate the sum ([U]) by the Cauchy identity 

oo / 

sx{x)sx{y) = n (1 - XiVj)'^ = exp I ^ Hkh 

A i,j=l \k=l 

where 

oo 

k 



k=l k=l 

Letting x = y = qP amounts to setting 

k/2 



k 



k ^ k(\ — q^ 
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Consequently, 



°° ^mk+nk+k \ 

^^'—) 

k=l m,n=0 / 

= exp(- log(l-g'"+"+^) ) 

\ m,n=0 / 

OO 

m,n=0 

Grouping the terms in the last infinite product with respect to the value of 
I = m + n + 1, we find that the partition function becomes the so called 
MacMahon function: 

OO 

Z = ll(l-qV- (8) 

1=1 



2.4 Models with external potentials 

The foregoing melting crystal model can be deformed by external potentials 
that depend on the main diagonal slice A = 7r(0). 

For example, we can insert the term Q'^' (Q > 0), namely, the potential 
|A| with coupling constant logQ- The partition function 

can be calculated in much the same way and becomes the deformed MacMa- 
hon function 

OO 

Z{Q) = l[{l-QqT'- (9) 

1=1 

We can further consider the potentials 

OO OO 

1=1 i=l 

which depend on an integer parameter p as well. These potentials originate 
in 5D supersymmetric U{1) Yang-Mills theory [2], [3]. Note that this definition 
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is rather heuristic; the infinite sums on the right hand side are divergent for 
the parameter q in the range < g < 1. A true definition is obtained by by 
pairing the term e'^^'^'"*"'"^ in the first sum with the term e^~*"^^ as 

$,(A,p) = y^^q^^P+^^-^+^) - + q^—^. (10) 

i=l 

The second term on the right hand side amounts to the terms that cannot 
be paired. Since \i = i for all but a finite number of i's, the sum on the right 
hand side is a finite sum. Thus $j(A,p) turns out to be well defined. 

Our previous work deals with the model deformed by the linear com- 
bination 

oo 

fc=i 

of these potentials with an infinite number of coupling constants t = (ti, ^2, • • •)• 
Its partition function reads 

Z(t,p,Q) = ^gklQk(0)l+p(p+l)/2g*{t.^(0).p)_ 
7r 

Note that the potential |A| of the weight Q''^' is also modified to |A|-|-p(p-|-l)/2. 
Unlike Z and Z[Q), this partition function cannot be calculated in a closed 
form. We could, however, show that it coincides, up to a simple factor, with 
a tau function of the ID Toda hierarchy. 

2.5 Models with two g-parameters 

We now turn to the models with two g-parameters [IT]. These models are 
obtained by substituting the Schur function factors as 

Sxiq'Y — s,iq',)sM)- 

The new g-parameters gi, g2 are assumed to be in the range < gi < 1 and 
< g2 < 1. This substitution amounts to modifying the g-weights q^, q^ of 
semi-standard tableaux T, T' on the main diagonal slice as 

rj-f rp rpt rpl 

q — 'Qi, q — ' ^2 ■ 
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The previous partition functions Z, Z{Q) and Z{t,Q,p) are thereby re- 
placed by 

Z{qi,q2) = ^Sx{qi)sx{q2), 

X 
A 

Z{t,Q,p;quq2) = ^AlgO^AfeOQ^'e*^*'"'^^- 

A 

Remember that the potentials ^{t,X,p) contain the third g-parameter q as 
well. It is these partition functions that we shall consider in detail. As 
regards the first two, we can apply the previous method to give an infinite 
product formula: 

oo 

ziquq2)= n ii - qr''\r''r\ 

'"r (11) 

ZiQ;quq2)= n (1-Q?r%2^'^')-^ 

m,n=0 

3 Fermionic formula of partition function 

We use a 2D complex free fermion system to reformulate the partition func- 
tions. Notations and conventions are the same as those in our previous work 

m- 

3.1 Complex fermions 

Let 'ipny'ipn ('^ ^ 2) denote the Fourier modes of the free fermion fields 

oo oo 

n=— oo n=— oo 

They satisfy the anti-commutation relations 

The Fock space splits into charge p subspaces J-'p {p G Z). JF^ has a 
normalized ground state (charge p vacuum) \p), which is characterized by 
the vacuum condition 

i^nlp) = for n > —p, i^nlp) = n > p + 1. 
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The dual Fock space JF*, too, splits into charge p subspaces T* (p G Z). The 
vacuum condition for the normalized ground state (p| of T* reads 

(pl'^n = for n < —p — 1, {p\4'* = for n < p. 

Tp and T* have excited states |A,p) and (A,p| labeled by partition A, which 
altogether form a basis of and JF^. If A is of length < n, namely, A = 
(Ai, A2, . . . , A„, 0, 0, . . .), these excited states can be obtained from \p) as 

{\,p\ = • ■ •^-(p-„)-l^(;+A„-n)+l ■ ■ ■^(p+Ai-l)+l- 

These states are mutually orthonormal: 

{X,p\/X,q) = Sxf,5pq. 



3.2 Partition functions of simplest model 

Let us introduce the special fermion bilinears 

00 00 00 

Jm= ■'^rn-n'ipn-^ ^0 = n:^_„<:, = ^ g'^^i^'-nC:, 

n=— 00 n=— 00 n=— 00 

where : : stands for normal ordering, namely. 



:^m<:=V'n^C-(0|^n^^:|0). 

Jm's are the Fourier modes of the U{1) current 



-m-l 

1 



and satisfy the commutation relations 

of the Heisenberg algebra. Lq is one of the basis {Lm}m=-oo of ^^e Virasoro 
algebra. H^s are the same "Hamiltonians" as used in the case of gi = g2 = 1 
[To] . HkS and Lq commutate with each other, and have the excited states 
{X,p\ and \X,p) as joint eigenstates, namely, 

{X,p\H, = (A,p|$fe(A,p)|, (A,p|Lo = (A,p| (\X\ + P^ll^) (13) 
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and 



Hk\\p) = ^k{X,p)\X,p), Lo|A,p)=(|A| + ^^^^)|A,p). (14) 



In particular, the potentials $fc(A,p) and their linear combination X,p) 
show up here as the eigenvalues of Hj. and their linear combinations 

oo 
k=l 

We further introduce the exponential operators [T] 

which belong to the Clifford group GL(oo). These operators can be factorized 
as 

oo 

G± = n V±ir^'^'), (15) 

m=0 

where V±{z) denote the familiar vertex operators 

V±{z) = exp X^T^^'^ ■ 



k 

These vertex operators act on (A,p| and \X,p) as 

{X,p\V4z) = V.{z)\X,p) =Y,z\'^\-\%,p). (16) 

One can thereby deduce [T] that the action of G± on the ground states {p\ 
and \p) yields a linear combination of {X,p\ and \X,p) with coefficients s\{q'^): 

A TishapcA A 

G-\P) = Y1 E g^'|A,p) = 5^s,(gO|A,p). 

A r':shapeA A 

G± thus play the role of "transfer (or transition) matrices" in lattice models. 
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Since (A,p| and \X,p) are orthonormal, the inner product of {p\G+ and 
G_\p) becomes a sum of sx{q^)^, 

A, /I A 

which is independent of p. We thus obtain the fermionic formula 

Z=(0|G+G_|0) (18) 

of the partition function Z. The MacMahon function formula can be also 
derived from from this formula and the commutation relation 

(oo k \ 

J^kd^^j (19) 

of G±. 

3.3 Partition functions of deformed models 

A fermionic formula of Z{Q) and Z{t,Q,p) can be obtained by inserting a 
combination of Lq and H^s with coupling constants into ( |T8l) . As regards 
Z{Q), the operator to be inserted is Q^°, which acts on (A, 0| and |A, 0) as 

(A,o|Q^« = (A,o|gi^i g^«|A,o) = gi^i|A,o). 

{0\G-^Q^°G-\0) can be thereby expanded as 



(0|G+g^«G_|0) = 5^(0|G+g^«|A,0)(A,0|G_|0) = 5^s,(g'')2gl' 



A A 



This is nothing but Z{Q). Thus we obtain the fermionic formula 

Z{Q) = (0|G+g^°G_|0). (20) 
In the same way, we can derive the fermionic formula 

Z(t,g,p) = (p|G+g^°e^(*)G_b) (21) 

for Z{t, Q,p). 

These fermionic formulae can be readily generahzed to Z(gi, 52), Z{Q]qi,q2) 
and Z{t, Q,p; qi, ^2)- Let us introduce the operators 

/ 00 k/2 

G+{qi) = exp ^ Jfc 



^ k{l - ql] 



00 k/2 

G_(g2) = exp ( ^ 
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{ k{l - ql 



in place of G±. As shown in ffTTj) . these operators generate a linear com- 
bination of {X,p\ and \X,p) with coefficients Sxiql) ^xil^.)- The rest of 
calculations in the previous case applies to this case without modification. 
Thus we obtain the following fermionic formulae of Z{qi, ^2), Z{Q; qi, ^2) and 
Z{t,Q,p; gi,g2): 

Z{quq2) = {0\G4q^)G^{q2m, 
Z{Q;qi,q2) = (0|G+(gi)g^«G_(g2)|0), (22) 
Zit,Q,p;q,,q2) = (p|G+(gi)Q^«e^(*)G_(g2)|p). 



4 Relation to Toda hierarchy 

We now turn to the issues of integrable structures. It is shown in our previous 
paper [lOj that the partition function Z{t,Q,p) coincides, up to a simple 
factor, with a special tau function of the ID Toda hierarchy. The goal of this 
section is to generalize this result to the case where qi and q2 are related to 
q as 

gi = g^/^\ 92 = g^/^^ (23) 

for some positive integers A^i and A^2- We assume this condition throughout 
this section. 



4.1 Intertwining relations in quantum torus algebra 

Let us recall yet another set of fermion bilinears 

00 00 

n=— oo n=— oo 

from our previous work [lOj. Wo is one of the basis {Wm}m=-oo of ^^e so 
called VF3 subalgebra in the W^o algebra of complex fermions. Vm ^'^s satisfy 
the commutation relations 

k+l 

[Vi'\ Ki')] = (g('™-^")/^ - g(^"-'™)/^)(Ki'i'^ - 5.+„,oY^) (24) 

of the quantum torus algebra]^ Actually, we need just a half of this algebra 
(so to speak, a quantum cylinder algebra) spanned by Vm '' with k > 1 and 
m G Z. Jm and Hk are part of this basis: 

Jm = Vi'\ H, = vt\ (25) 

^Substantially the same realization of this algebra is considered in different contexts 
by Gao [25 and Okounkov and Pandharipande [3 IB] ■ 
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This Lie algebra has an inner symmetry ("shift symmetry") [TU]. A conse- 
quence of this symmetry is the following intertwining relations among J^'s 
and HkS, which play a fundamental role in identifying the integrable struc- 
ture of Z(t, Q,p). 

Lemma 1. Jm's and H^'s satisfy the intertwining relations 



(26) 



We now use these intertwining relations replacing q qi, q2- f !23l) implies 
that Hk may be thought of as elements of the quantum torus algebras with 
g-parameters q = qi and q = q2, namely, 

Hk = V,^''^'\q,) = vi''^'\q2), (27) 

where Vm\qi) and Vm\q2) denote the counterparts of Vm^ for q = qi and 
q = q2- Therefore, applying (!26l) to the cases where q = qi and g = g2, we 
obtain the relations 

qY'/'G4q^)G^iq^)Hk= (^i-l f^'jN.k + qT"^^^ (28) 

interchanging JNik^ and H^'s and 

HkG4q2)G^iq2)qY'^' = G^{q2)G^{q2)q'^'" {{-I)''-' J^^.k + (29) 

interchanging J^n^^s and HkS. Note that the c-number terms in the paren- 
theses have been rewritten as 

Nik N2k k 

(30) 



1 - g^i'^ 1 - q^^f" 1-q 

4.2 Partition function as tau function 

Armed with these intertwining relations, we can generalize our previous 
result pro] to Z{t,Q,p;qi,q2). Let us recall that a general tau function 
t{T,T,p) of the Toda hierarchy depends on two sets of continuous vari- 
ables T = (Ti, T2, . . .) and T = (Ti, T2, . . .) and has a fermionic formula of 
the form 



t{T,T,p) = (p|exp Ij^TkJk gexp l-J^TkJ-k \p), (31) 



,fc=l / \ k=l 
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where g is an element of GL{oo). We now show that Z{t,Q,p;qi,q2) coin- 
cides, up to a simple factor, with such a tau function. 

Theorem 1. Z{t,Q,p;qi,q2) can be rewritten in two different forms as 

qH, 



Z{t, Q.p, qu q2) = [q^q^r^^^^'^^'^^'^'' exp ^ 

(oo 
Y^[-lf^%JN., ] g\p) (32) 
k=l 

and 

(oo ^ 
fe=l ^ 

/ oo 

X (pl^^exp Y,i-^f''tkJ-N,k ) b), (33) 



.fc=i 



where g is an element of GL{oo) of the form 

g = qY''^'G4q,)G+iq^)Q'^^G4q2)G^iq2)qY'^'. (34) 



Proof. To derive flHzj) . let us think of the right hand side of the fermionic 
formula ( l22l) as the inner product of (p|G+(gi) G JF* and G-{q2)\p) G in 
which Q^°e^^*^ is inserted. Also remember that the order of and e^*^*-* 
is immaterial, because Lq and H{t) commute with each other. Since G'-(q'i) 
and qY"^"^ act on {p\ almost trivially as 

(plG^qi) = {Pl (p|g^°/^ = gf^-^^)(^^+^)/^^(p|, 
we can rewrite (p|G'+(gi) in ( l22l) as 

(pIG^gi) = gr^^^^^^^^^^^^/"(p|gr/^G_(gi)G^(gi). 

In the same way, 

G-i^M = g2-^(^^^^^^^^^^/"G_(g2)G^(^2)gr"/^b). 

Thus Z{t, Q,p; qi, q2) can be cast into such a form as 

Z(t,g,p;g„g2) = (gig2)-^(^+^)(^^-^^)/^^ 

X (p|gr''^'G_(gi)G+(gi)e^WQ^°G_(g2)G+(g2)gr"/'b). 
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Since the intertwining relation implies the identity 



we can move e^^*^ to the right of {p\ and obtain the first formula (|32|) . The 
second formula flHHj) can be derived in much the same way. □ 

Thus Z{t, Q, p; qi, 52) turns out to coincide, up to a simple factor, with the 
tau function t{T,T,p) determined by the GL(oo) element g of flMl) . Note 
that the time variables T, T of the full Toda hierarchy are now restricted to 
a subspace. In particular, unlike the case where iVi = A''2 = 1, not all of these 
variables join the game. Namely, it is only Tat^/j and T/Vafc, k = 1,2, . . ., that 
correspond to the coupling constants t of this generalized melting crystal 
model. 



4.3 Bigraded Toda hierarchy 

Since the same partition function Z{t, Q,p; qi, ^2) is expressed in two appar- 
ently different forms as flH^ and flHHl) . we find that the identity 



(p|exp \^(,~lf^%JM,,j g\p) = {p\gexp J_^^, j \p) (35) 

holds. This is a manifestation of the following more fundamental fact. 

Theorem 2. Jk 's and the GL{oo) element g of [34\ ) satisfy the intertwining 
relations 

{-lf-KjN,kg = g{-lf''J-N,k (36) 

fork = l,2,.... 



Proof. Using the intertwining relations fl28p and fl29l) . we can derive (1361) as 
follows: 

{-l)''''jN,kg = {-lf^'JN,kqT''^'G4q^)G^iq^)Q'^^>G4q2)G^iq2)qY"^' 

= gr°/'G_(gOG+(gi) (h, - Q^^^G.iq,)GM'h''^' 

= qY'"G^{q,)G+{q,)Q'^^ [h, - G.{q,)G^{q,)qf^" 



q7'''G4q,)G4q,)Q^^G4q,)G^iq2)qr'\-ir''J-N,k 
16 



gi-lf-'j- 



□ 

By these intertwining relations, we can freely move Jn^u and J-N2k to the 
far side of g. This implies that t{T,T,p) depends on T/v^fc and T/V2fc only 
through the linear combination {—l)'^'^^T]^-^^k — {—l)'^^^Ti^2k. In other words, 
t{T,T,p) satisfies the constraints 

hir'^ + {-ir'^ = (37) 

for k = 1,2,.... Apart from the presence of the signature factors (— l)^i^ 
and (— 1)^2^^, these constraints are the same as those that characterize the 
bigraded Toda hierarchy of type (A^i, A^2) [13 as a reduction of the full Toda 
hierarchy. 

In the language of the Lax operators [I3j 

L = e^" +ui + U2e'^'' H , L = uoe^" + uie^^" H 

of the Toda hierarchy, the reduction to the (slightly modified) bigraded Toda 
hierarchy can be characterized by the constraint 

(-L)^i = (-1)-^^ (38) 

Let C denote the difference operator defined by both hand sides of this con- 
straint. This reduced Lax operator is a difference operator of the form 

C = (-e^.)N, + b^i-e'^f^-' + ■■■ + 6^^^.^^(-e^^)-^^ (39) 

and satisfies the Lax equations 

§=|ft.£|. = (40) 

where and B^. are given by 



Bk = {L'')>o, Bk = {L ^] 



<05 



( )>o and ( )<o standing for the projection onto the part of nonnegative 
and negative powers of e^^ . More precisely, since the ( )>o and the ( )<o 
parts of £ = {—L)^^ = {—L)~^^ are given by 

iC)>o = i-lf'B^,, (/:)<o = (-l)^^5;v2, 
we can express £ itself as 

C = {-lf^BN, + i~-lf^BN,. (41) 
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The powers £ , k = 1,2, . . ., of C can be hkewise expressed as 
This readily implies that the stationary equations 

OlNik OlN2k 

hold for C These stationary equations are counterparts of (!37|) in the Lax 
formalism. 



5 Relation to g-difFerence Toda equation 

In this section, we consider another integrable structure, which is hidden in 
Q-dependence of the partition function. We no longer have to assume ( l23ll . 
namely, qi, q2 and q are independent parameters throughout this section. 

5.1 Variables for g-difFerence analogue 

It is known pOl EU |22] that tau functions of the Toda hierarchy can be 
converted to tau functions of a g-difference analogue of the 2D Toda equation 
by changing variables as |fl 

k k 

k{l - ql) k[l - gf) 

Let a{x,y,p) denote such a transformed tau function, namely, 

(^{x,y,p) = T{[x]q,,-[y]q2,p), (44) 

where t(T, T,p) is a tau function of the Toda hierarchy and [x]q denotes the 
g- difference analogue 

9 k 



l-q 2{l-q^y'' k{l-qk)' 



of the notation 



^2 

\X\ I X, 2 ) • • • ; 



■^This definition is slightly modified from the common one in the literature [201 [^[^ . 
We shall make a remark on this issue later on. 
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that play a fundamental role in the study of KP and Toda hierarchies. 

By this change of variables fHHj) . the fermionic formula (IHTj) of tau func- 
tions of the Toda hierarchy turns into the formula 



a{x,y,p) = (plexp (^g ^^^^ J. j ^exp (^g j \p) (45) 

of a{x,y,p). In particular, when x,y and g are specialized as 

x = g^/^ y = gr^ 9 = Q''e^'% (46) 

a{x, y,p) coincides with Z{t, Q,p; qi, 52)- The goal of this section is to derive 
a g-difference equation that Z(t, Q,p; qi, ^2) satisfies with respect to Q. 

5.2 g-difference 2D Toda equation for a{x,y,p) 

As a preliminary step towards a g-difference equation for Z{t, Q,p; qi, ^2), we 
now show that the transformed tau function a{x,y,t) satisfies a g-difference 
analogue of the 2D Toda equation (in a bilinear form) 

d^T{x,y,p) ^ y p) - "^^(^'^'P) ^'^(^^y^P) 
dxdy ' ' dx dy 

= T{x,y,p + l)T{x,y,p-l). (47) 

Now that this is the lowest equation satisfied by the tau function of the Toda 
hierarchy with respect to x = Ti and y = —Ti. 

Lemma 2. For any tau function of the Toda hierarchy, the function cr{x, y,p) 
defined by [44^ satisfies the q-difference 2D Toda equation 

(^{qix, q2y, p)(^{x, y, p) - (y{x, q2y, p)(T{qiX, y, p) 

= xya{x, y,p+ l)a(gia;, q2y, p-l). (48) 

Proof This is a consequence of the difference analogue 

r(T-[x],T,p)T{T,T-[y],p)-T{T,T,p)r{T-[x],T-[y],p) 

= xyriT, T~[ylp+ l)r(T - [x], T, p - 1) (49) 

of the 2D Toda equation ( 1471) . This equation is one of "Fay-type identities" 
that hold for any tau function of the Toda hierarchy. We shift T as 
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T ^ T + [y] and substitute T = [x]g^ and T = — [y]g2 in this equation. The 
outcome is the equation 

T{[x]g, - [x],-[y]g^ + [y],p)T{[x]g„-[y]g^,p) 

-r{[x]g^,-[y]g,^ + [y],p)T{[x]g^ - [x],-[y]g^,p) 
= xyT{[x]g„ - [y]g„p + l)T{[x]g^ - [x],-[y]g, + [y],p-l), 

which we can further rewrite as 

- T{[x]g,, -[q2y]q2,PM[qix]g„ -[y]g2,p) 

= XyT{[x]g^, ~[y]q2^P+ '^)'r{[qiX]q,,-[q2y]q2,P- 1) 

by the g-shift property 

[qx]g = [x]g - [x] (50) 
of [x]g. The last equation is nothing but ( l48l) . □ 
A few remarks are in order. 

1. One can rewrite the g-difference equation ( HHl) as 

Dqi,xDg^^ya{x, y, p) ■ a{x, y, p) - y, p) ■ Dg^^ya{x, y, p) 

= a{x,y,p+l)a{qi,q2y,p- 1), (51) 

where -Dgi,a: and Dg^ y stand for the g-difference operators 

(^ix,y,p) - a{qix,y,p) 



Dqi,x(^{x,y,p) 
Dq2,y(^{x,y,p) 



(l-gi)x 
(^ix,y,p) - a{x,q2y,p) 



(1 - q2)y 

As gi, g2 0, this equation turns into the 2D Toda equation (H7I) . 
2. One can modify fHHj) as 



k h 

X - y 

^' = W^y ^' = w^y ^^^^ 

The transformed tau function 

p{x,y,p) = T{[x]g„[y]q,,p) (53) 
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has a fermionic formula of the form 

p{x,y,p) = (plexp (^g^^^J.j ge^v g ^(^^^-.j b)(54) 

and satisfies the equation 

p{qix, y, p)p{x, q2y, p) - p{x, y, p)p{qiX, gaZ/, p) 

= xyp{x,q2y,p+l)p{qix,y,p- 1). (55) 

This is another version of the g-difference Toda equation, which is 
rather common in the hterature [201 EH ES] . 

3. One can consider the functions 

oo 

^{x,y,p) = a{x,y,p) JJ (I - xyq'^'q^), 

(56) 



oo 

m „n\ — 1 



p{x,y,p) = p{x,y,p) Yl{i-xyq'^q2 

m,n=0 

in place of a{x, y, p) and p(x, p). In the usual formulation of the Toda 
hierarchy, these modified tau functions amounts to 

f (T, T, p) = t{T, T, p) (|- E ^^'^^'^j ' (5'^) 

which is suited for the construction of the so called Wronskian solutions 
of the Toda hierarchy. a{x,y,p) and p{x,y,p) satisfy the g-difference 
equations 

^{qix, q2y, p)&ix, y, p) - (1 - xy)cr{x, gaZ/, p)^{qix, y, p) 

= xya{x,y,p + l)a{qix,q2y,p-l) (58) 

and 

p(gix, y, p)p{x, gay, p) - (1 - xy)p{x, y, p)p{qix, gsy, p) 

= xyp{x,q2y,p+l)piqix,y,p- 1). (59) 

The equation for p{x, y, p) coincides with the g-difference Toda equation 
of Kajiwara and Satsuma 



21 



5.3 g-difference equation for Q, j9; gi, ^2) 

As already mentioned in the beginning of this section, Z{t,Q,p;qi,q2) is 
equal to a special value of the function a{x,y,p) defined by the fermionic 
formula ( H5|) with g = Q^°e^^*\ namely, 

Z{t,Q,p;qi,q2) = a{ql^^ ,ql^^ ,p). (60) 

To derive a g-difference equation with respect to Q, it is more convenient to 
consider this relation in a slightly more general form as follows. 

Lemma 3. This special solution a{x,y,p) of the q-difference Toda equation 
is related to the partition function Z{t, Q,p; qi, ^2) 

(^{x,y,p) = {qi^^^q2^^'^xyyP'-P^'^^^^Z{t,q^^^^q2^^^xyQ,p;qi,q2). (61) 
Proof. Since Lq and J^s satisfy the well known commutation relations 

we can rewrite the exponential operators in (H5ll as 



(°° k \ / 00 1/2 \ 

Consequently, 

/ 00 1/2 \ 
a{x,y,p) = (j9|(gr'/'x)-^°exp ( j (g-V2^)^ogLo^//W 

X (g2-^^)"«exp ( E fc(ri:^^-0 (92"'^)-"lp). 



Since 



(p|(g-'/\)-^« = (gr'^'x)-*'(p+^)/2(p|, 
(g2"%)-^1j9) = (g2"^/'l/)-^(^+^)/'b) 



and 



(gr^/^x)^"g^"e^(*)(g-^/^y)^" = (gr'/^g2"'/^a;yQ)^0e^W, 
this expression of a{x, y,p) boils down to flUTl) . □ 

We can now derive a g-difference equation for Z(t, Q, ; gi, ^2) from the 
g-difference 2D Toda equation (HHl) as follows. 
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Theorem 3. Z{t,Q,p]qi,q2) satisfies the q-difference equation 
Z{t, qiq2Q; qi, q2)Z{t, Q, ; gi, gs) - Z{t, qiQ, ; gi, q2)Z{t, gsQ, ; Qi, qi) 

= (gig2r'/'^(t,g,P+l;gi,g2)^(t,gig2Q,p-l;gi,g2). (62) 

Proof. Plugging a{x,y,p) of ( 16T|) into ( 148|) yields the equation 

Z{t, q^^ql^^xyQ, ; qi, q2)Z{t, q^^''^q2^''^xyQ; qi, ^2) 

- Z{t, ql^'^q2^^^xyQ, ; qi, q2)Z{t, qi^^'^ql^^xyQ; qi, 52) 
= {(li(l2T^^^^Z{t, qi^^^q2^^^xyQ,p+ 1; gi, g2)^(*, ql^'^ql^^xyQ,p-l; qi, ^2)- 

Upon rescaling Q as Q x~'^y~^qy^ql^'^Q, this equation turns into flU^ . □ 

It will be instructive to compare this result with the reduction process 
from the 2D Toda equation (H7|) to the ID Toda equation 

- (^)' = rit,,^l)r(t,p - 1). (63) 

(l63l) is obtained from ( 1471) by assuming the condition that 

^{x,y,p) = T{x + y,p), 

namely, T{x,y,p) be a function of t = x + y and p. In the same sense, if 
cr(x, y,p) is a function of t = xy and p, namely, 

(^{x,y,p) = (r{xy,p), 

then the g-difference 2D Toda equation ( l48i) reduces to the g-difference ana- 
logue 

o-(?i?2t, J3)o-(t, P) - cr(git, p)cr(g2t, p) = tcr(t, J9 + l)cr(gig2t, p - 1) (64) 

of the ID Toda equation ( l63i) . The g-difference equation ( l62l) stems from 
substantially the same idea. The apparent discrepancy of (!62|) and (IMI) is 
due to the prefactor [q^^^"^ q^^^"^ xy)^^^^^^^^ on the right hand side of fl6Tl) . 

6 Conclusion 

We have found two kinds of integrable structures hidden in the partition 
function Z{t,Q,p]qi,q2) of the generalized melting crystal model. The first 
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integrable structure is the bigraded Toda hierarchy that determines the t- 
dependence of the partition function. This integrable structure emerges when 
gi,g2 and q satisfy the algebraic relations ( !23l) (Theorems 1 and 2). This is 
a natural generalization of the main result of our previous work [10]. The 
second integrable structure is a g-difference analogue of the ID Toda equation 
(Theorem 3). The role of time variable therein is played by Q, and this 
equation holds without any condition on qi, q2 and q 

It is easy to see that these results still hold if a GL{oo) element of the 
form 

(oo 
^ Cn'.lp-ni'n- 
n=— oo / 

is inserted in front of Q^" (though the ordering of h, and e^*^*^ can be 
changed arbitrarily). Of particular interest is the case where 

h = e^^" {P is a constant). 

This amounts to introducing a new potential of the form @ 



J](p+A.-. + l)^-5^(-. + l)^ 

1=1 i=l 

= f:((p+A.-.+ l)^-(p-. + lrt + £<P±IK^ (65) 

i=l 

with coupling constant (3. The partition function 

Z{t,(3,Q,p;q„q2) = (p|G+(gi)e^^«g^V^WG_(g2)b). (66) 

thus obtained can be found in 5D supersymmetric Yang-Mills theory [21 [3] , 
and also related to topological strings on a class of toric Calabi-Yau threefolds 
and Huwitz numbers of the Riemann sphere [30] . 
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